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ABSTOACT  Protective  antigen  is  essential  for 
the  pathology  of  Bacillus  anthracis  and  is  the  pro¬ 
posed  immunogen  for  an  improved  human  anthrax 
vaccine.  Known  since  discovery  to  comprise  differ¬ 
entially  charged  isoforms,  the  cause  of  heterogene¬ 
ity  has  eluded  specific  structural  definition  until 
now.  Recombinant  protective  antigen  (rPA)  con¬ 
tains  similar  isoforms  that  appear  early  in  fermenta- 
tion  and  are  mostly  removed  through  purification. 
By  liquid  chromatography-tandem  mass  spectrome¬ 
try  sequencing  of  the  entire  protein  and  inspection 
of  spectral  data  for  amino  acid  modifications,  phar¬ 
maceutical  rPA  contained  measurable  deamidation 
at  seven  of  its  68  asparagine  residues.  A  direct  asso¬ 
ciation  between  isoform  complexity  and  percent 
deamidation  was  observed  such  that  each  decreased 

-lth^ty  and  lncreased  with  Protein  aging.  Posi¬ 
tion  N537  consistently  showed  the  highest  level  of 
modification,  although  its  predicted  rate  of  deami¬ 
dation  ranked  10th  by  theoretical  calculation,  and 
other  asparagines  of  higher  predicted  rates  were 
observed  to  be  unmodified.  rPA  with  more  isoforms 
and  greater  deamidation  displayed  lower  activities 
for  furin  cleavage,  heptamerization,  and  holotoxin 
formation.  Lethal  factor-mediated  macrophage 
toxicity  correlated  inversely  with  deamidation  at  re¬ 
sidues  N466  and  N408.  The  described  method  meas- 
ures  deamidation  without  employing  theoretical 
isotopic  distributions,  comparison  between  differ¬ 
entially  treated  samples  or  computational  predic¬ 
tions  of  reactivity  rates,  and  is  broadly  applicable  to 
the  characterization  of  other  deamidated  proteins. 
Proteins  2007 ;68:458-479.  ©  2007  Wiley-Liss,  Inc. 
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INTRODUCTION 

The  protective  antigen  (PA)  of  Bacillus  anthracis  plays 
an  essentia]  role  in  anthrax  pathology,  and  is  under 
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mic  transporter  component  for  two  catalytic  moieties,  the  documentation  of  methods  and  results  that  provide 

edema  factor  (EF),  which  is  a  calmodulin-sensitive  ade-  increased  understanding  of  its  physicochemical  and  mo- 

nylate  cyclase,  and  lethal  factor  (LF),  a  Zn2+-dependent  lecular  structure.29’30  With  regard  to  biological  activity 

endoproteinase  of  mitogen-activated  protein  kinase  ki-  in  vivo ,  prior  investigation  has  shown  that  the  two  major 

nases.5-8  The  deregulation  of  cellular  physiology  by  EF  isoforms  and  principal  components  of  pharmaceutical 

and  LF  leads  to  tissue  swelling,  necrosis,  shock,  and  grade  rPA  elicit  equivalent  specific  immunity  and  protec- 

rapid  death  if  untreated,  all  hallmark  symptoms  of  an-  tion  against  lethal  challenge  by  B.  anthracis  spores, 

thrax  9  PA  protein  (83  kDa)  binds  to  cell-surface  anthrax  The  analytical  method  and  findings  reported  here  reveal 

toxin  receptors  (ATR)10  and  is  cleaved  by  the  furin  fam-  multiple  site  deamidation  as  the  chemical  source  of 

ily  of  proteases11  to  create  a  63  kDa  active  form  and  a  micro-heterogeneity  that  was  long  observed  as  nega- 

20  kDa  fragment  of  undefined  function.12  Cleaved  PAea  tively  charged  isoforms  in  native  PA  protein, 

assembles  into  a  heptameric  prepore  (PA«3)7  that  binds  Purified  proteins  are  particularly  susceptible  to  degra- 
up  to  three  molecules  of  EF  or  LF.1314  This  complex  is  dation  through  the  nonenzymatic  post-translational  mod- 
endocytosed  and  trafficked  to  an  acidic  compartment  ifications  of  deamidation,  oxidation,  and  aggregation, 
wherein  a  large  loop  from  domain  2  in  each  PA  subunit  Because  any  of  these  modifications  may  affect  the  corn- 

inserts  into  the  membrane  and  forms  a  14-strand  [1-bar-  position,  quality,  consistency,  or  stability  of  the  final 

rel.15~18  A  model  of  translocation19'20  proposes  that  the  form,  degradation  activities  introduce  the  need  for  spe- 

EF  (89  kDa)  and  LF  (90  kDa)  pass  through  the  channel  cial  controls  during  the  manufacture,  storage,  and  use  of 

of  this  pore  This  model  invokes  unfolding  of  EF  and  LF  purified  proteins  intended  for  scientific  study  or  com- 

for  transit,  with  spontaneous  refolding  in  the  target  cell’s  merce.  Deamidation  is  generally  the  most  prevalent  of 

cytosol,  since  the  X-ray  structures  of  EF  and  LF  have  these  spontaneous  changes  to  protein  structure.  The 

minimum  measured  diameter  of  -45  A,  while  the  PA  deamidation  reaction,  which  is  the  loss  of  ammonia  from 

pore  has  a  diameter  of  ~12  to  -17  A.18  Alternatively,  the  amide  side  group  of  asparagine  (Asn)  or  glutamine 

full  endoproteinase  activity  of  LF  was  demonstrated  to  (Gin),  proceeds  through  a  cyclic  imide  intermediate  (see 

occur  without  dissociation  from  the  membrane-bound  PA  Fig.  1)  whose  rate  of  formation  and  hydrolysis  is  strongly 

complex.21  Although  the  precise  mechanism  is  still  affected  by  the  adjoining  distal  amino  acid,  local  spatial 

uncertain,  it  is  clear  that  factor  binding  is  essential  and  structure,  and  solvation  chemistry,  including  pH  and  sol- 

the  net  effect  of  PA  activity  is  to  expose  targeted  host  ute  composition. 33-39  Deamidation  introduces  an  ioniz- 

cells  to  the  catalytic  cytotoxic  factors.  By  itself,  however,  able  carboxylate  group,  which  imparts  a  negative  charge 

PA  is  nontoxic  and  elicits  protective  immunity  against  B.  at  physiological  pH  and  increases  the  mass  of  the  pro- 

anthracis,  as  demonstrated  by  numerous  animal  stud-  tein  by  1  Da.  Asparagine  residues  were  first  observed  to 

ies  2,22,28  ’Because  of  its  pivotal  role  in  pathogenesis  and  be  -40  times  more  prone  to  deamidation  than  Gin  resi- 

its  immunoprotective  activity,  PA  has  accordingly  been  dues  within  proteins.34  More  recently,  extensive  studies 

an  unknown  or  intended  principal  constituent  in  an-  of  model  synthetic  peptides  indicate  that  Asn  residues 

thrax  vaccines  since  the  first  demonstration  of  protective  deamidated  on  average  about  100  times  faster  than  do 

preparations  125  years  ago.1'22  Gin  residues  in  peptides.  Consequently,  most  expen- 

The  modern  search  for  a  safe  and  effective  human  an-  mental  observations  concern  Asn  deamidation  and  com- 

thrax  vaccine  has  led  to  the  development  of  purified  monly  measure  one  of  its  two  reaction  products, 

recombinant  PA  (rPA),3'23’24  considered  the  top  candidate  Although,  Asn  deamidation  yields  normal  L-aspartate 

for  the  first  subunit  vaccine  to  replace  existing  anthrax  (Asp),  the  more  abundant  product  is  isoaspartate 

vaccines,  which  are  comprised  of  alum-precipitated  B.  (isoAsp),  a  p-linked  isomer  which  incorporates  the  side 

anthracis  culture  supernatant  preserved  with  formalin.  chain  Cp  into  the  peptide  backbone  (see  Fig.  1). 

The  currently  licensed  human  anthrax  vaccines  (Bio-  Because  isoAsp  also  lengthens  the  polypeptide  chain  by 

Thrax®,  BioPort,  Lansing  MI;  and  Anthrax  Vaccine  Pre-  one  methylene  group  (CH2),  it  is  more  deleterious  to  pre¬ 
cipitated,  Porton  Down,  UK)  are  based  on  50-year-old  tein  structure  and  function, 32'35'39'41  and  consequently 

technology  are  reactogenic  and  require  multiple  doses  to  has  attracted  more  attention  in  the  characterization  of 

impart  full  immunity  with  annual  boosting  for  mainte-  protein  stability.31'42’43  This  aberrant  p-isopeptidyl  bond 

nance.1-3  In  addition,  they  contain  many  cellular  compo-  is  corrected  in  vivo  by  protein  L-isoaspartyl 

nents  besides  PA  of  unspecified  influence  on  product  O-methyltransferase  (PIMT,  EC2.1.1.77),  an  enzyme  first 

effectiveness  and  safety.25  Animal  studies  have  demon-  discovered  in  brain  tissue,  but  which  now  appears  to  be 

strated  that  purified  rPA  elicits  equivalent  immunoge-  widely  distributed  and  critical  for  alleviating  this  form  of 

nicity  and  protective  immunity  compared  with  these  protein  aging.36’41’44  A  growing  list  of  maladies  and  cel- 

licensed  vaccines.2'3’26-28  As  an  active  pharmaceutical  in-  lular  dysfunctions  are  associated  with  PIMT  deficiency, 

gredient  intended  for  eventual  use  in  humans,  rPA  increased  isoAsp  content,  or  Asn/Gln  deamidation, 

requires  thorough  characterization  in  many  aspects  of  including  cataracts,45’4  Alzheimer  s  disease,  spina 

product  safety  and  effectiveness,  including  detailed  bifida,49  autoimmunity,50  and  prion-related  encephalopa- 

descriptions  of  biological  activity  and  chemical  structure.  thies.51'52  Spontaneous  deamidation  has  also  been  associ- 

When  heterogeneity  is  evident  in  a  proposed  biological  ated  with  normal  processes,  such  as  Bcl-XL-regulated 

product,  as  is  the  case  for  rPA,  FDA  regulation  specifies  apoptosis53  and  homeostasis  of  cytochrome  C.  Because 
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caused  by  partial  and  differential  deamidation  across 
seven  Asn  residues,  and  to  observe  direct  correlations 
between  their  deamidation,  isoform  complexity,  and 
reduced  biochemical  and  biological  activities  of  protective 

antigen. 

MATERIALS  AND  METHODS 
Biological  and  Chemical  Materials 

Four  grades  of  rPA  protein  were  used  in  this  study. 
Pharmaceutical  grade  rPA  was  purified  from  ASterne- 
1(dPA102)  CR4,  a  nonsporogenic  avirulent  expression 
strain  of  B.  anthracis ,68’69  by  contract  (Biopharmaceut- 
ical  Development  Program,  SAIC-Fredenck,  MD)  un  er 
good  manufacturing  practices  (GMP),  using  a  proprie¬ 
tary  procedure  derived  from  the  discovery-based 
method.70  This  GMP  grade  rPA  was  supplied  in  vials  as 
bulk  substance  at  1.18  mg/mL  and  94.6%  purity,  in 
50  m M  ammonium  acetate  at  pH  8.9.  Intentiona  y 
degraded  rPA  protein  of  the  same  GMP  lot  was  achieved 
through  exposure  to  multiple  cycles  of  freeze  thaw,  and 
hold  at  4°C.  Research  and  development  (R  &  D)  prepara¬ 
tions  of  rPA  were  purified  from  culture  supernatants  of 
the  same  B.  anthracis  expression  strain,  essentially  as 
published,70  except  pools  differing  in  isoform  content 
were  collected  as  separate  fractions  across  the  primary 
peak  of  the  second  step  in  the  sequential  process,  and 
held  frozen  at  -70"C  until  this  study.  Commercial  grade 
rPA  was  purchased  as  a  lyophilized  preparation  and 
resuspended  in  water  just  before  use  *  *h' 

manufacturer  (List  Biologicals,  Campbell,  CA).  LF  was 
prepared  as  described.4  Sequencing-grade  trypsin  was 
from  Promega  (Madison,  WI).  Chymotrypsin  and  V8 
(endoproteinase  Glu-C)  were  from  Sigma  (St  Louis, 
MO)  Synthetic  peptides  were  from  New  England  Pepti¬ 
des  (Gardner,  MA),  and  comprised  the  native 
(IAFGFNPNGNLQYQGK)  and  fully  deamidated  PIM1- 
corrected  (IAFGFNPDGNLQYQGK)  mimetics  of  the  PA 
tryptic  peptide  529-545.  HPLC-grade  acetomtn  e  with 
formic  acid  (LC-MS  CHROMASOLVfy  and  all  other 
chemicals  were  from  Sigma.  Sterile,  deionized  water  was 
polished  in  the  laboratory  by  a  Solutions  2000  water  pu¬ 
rification  system  (Aqua  Solutions,  Jasper,  GA). 


Capillary  Isoelectric  Focusing 

The  isoelectric  point  (pi)  of  each  measurable  rPA  pro¬ 
tein  isoform  was  obtained  by  capillary  isoelectric  focus¬ 
ing  (C1EF)  using  reagents  and  software  of  the  BioFocus 
instrument  as  described  by  the  vendor  (Bio-Rad  Hercu¬ 
les  CA).  Relative  mobility  (Rm)  of  proteins  used  to  con¬ 
struct  the  standard  curve  was  calculated  by  using  the 
formula:  Rm  =  (tv-tcW*-tc),  with  anafyte  migration 
time  denoted  as  “fp”  for  protein  standard,  te  for  cathode 
interface,  and  “<a"  for  anode  interface.  Theoretical  pi  val¬ 
ues  were  calculated  from  the  substitutions  of  the  PA 
amino  acid  sequence  (P13423),71  as  adjusted  m  sdico  to 
omit  the  29  amino  acid  leader  sequence  using  the  pro¬ 
tein  parameters  tool  of  ExPASy  (www.expasy.org). 


Polyacrylamide  Gel  Electrophoresis 

rPA  was  fractionated  on  native  PhastGels  (Amersham 
Biosciences,  Piscataway,  NJ)  as  described  previously 
employing  multiple  loadings  of  2  pg  per  lane  to  purify 
isoforms  for  LC-MS/MS  analysis.  Protein  isoforms  were 
also  fractionated  by  two-dimensional  polyacrylamide  gel 
electrophoresis  (2-DE)  according  to  instructions  of  the 
instrument  and  supplies  vendor  (BioRad,  Hercules,  CA), 
loading  12  pg  of  rPA  protein  in  sample  dehydration 
buffer  per  duplicate  IEF  strips  of  pH  ranges  8-10.  4-7, 
and  4. 7-5. 9.  For  second  dimension  separation,  the  deve  - 
oped  IPG  strips  were  equilibrated  as  prescribed  then 
loaded  onto  Bis-Tris  XT  Criterion  gels  (4-12%  IPG  well 
+1)  processing  six  gels  simultaneously  (including  repli¬ 
cates)  in  a  DodecaCell  apparatus  at  200  V  for  45  mm 

with  MES  running  buffer  system.  . 

For  imaging  and  protein  excision,  gels  were  washed  in 
deionized  water  and  stained  with  GelCode  blue  (Pierce) 
according  to  vendor  instructions.  Images  were  captured 
with  a  FlouS  Max2  (BioRad,  Hercules,  CA)  and  spot  pat¬ 
terns  analyzed  using  PDQuest  version  7.1.  Orthologous 
protein  isoform  spots  from  replicate  gels  of  each  separa¬ 
tion  were  selected  and  excised  using  the  Proteome  Works 
Spot  Plus  Cutter  (BioRad,  Hercules,  CA),  or  with  a  clean 
razor  blade.  All  gel  plugs  were  processed  promptly  tor 
in-gel  proteolysis  and  LC-MS/MS  analysis. 


Mass  Spectrometry 

Complete  amino  acid  sequencing  by  LC-Mb/MH 

The  complete  amino  acid  sequence  of  GMP  grade  rPA 
protein  was  determined  from  analysis  of  spectral  data 
collected  by  LC-MS/MS,  as  described  previously  but  with 
some  modifications.72  Protein  was  digested  with  trypsin, 
chymotrypsin,  Glu-C,  or  a  combination  of  these  pro¬ 
teases,  performed  as  recommended  by  the  vendor  For 
protein  contained  within  excised  gel  pieces  each  gel 
piece  was  diced  into  smaller  pieces  (~1  mm  ),  collected 
into  one  tube,  and  washed  3  X  10  min  in  acetonitrile:  25 
m M  ammonium  bicarbonate  (1:1),  dried  under  vacuum, 
then  resuspended  in  25  pL  of  25  m M  ammonium  bicar¬ 
bonate  pH  8.4  to  which  was  added  75  pL  (300  ng)  of 
trypsin  stock  (Promega),  prepared  previously  by  sus¬ 
pending  the  enzyme  in  25  m M  ammonium  bicarbonate 
at  4  pg/mL  and  stored  at  -20  C  as  recommended.  Frag¬ 
ments  were  recovered  in  gel  supernatants  and  combined 
with  two  consecutive  extractions  of  the  gel  pieces  with 
100  pL  of  50%  acetonitrile,  0.1%  formic  acid.  The  pooled 
extracts  were  dried  under  speed  vacuum  mid  reconsti¬ 
tuted  in  30  pL  of  5%  B  chromatography  buffer.  Samples 
were  transferred  to  1.1-mL  crimp-capped  vials  with  400- 
uL  volume  inserts  (PN  9301-1388,  Agilent  Technologies, 
Wilmington,  DE)  and  stored  at  4°C  for  analysis  by  LC- 
MS/MS  within  24  h.  For  direct  proteolysis  of  free  prot™ 
in  solution,  the  process  was  carried  out  as  recommended 
bv  protease  vendors  (18  h,  37°C)  using  11-15  pg  of  rPA 
protein  (enzyme/protein  ratio  of  1:15-20,  w/w).  The  reac¬ 
tion  was  stopped  with  4  pL  of  formic  acid.  Final  injected 
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amount  of  28  gL  contained  a  calculated  4.96-5.43  pmol 
digested  protein.  Importantly,  since  rPA  contains  no 
cysteines  protein  and  peptide  samples  were  not  dena 
tured,  reduced,  or  alkylated  in  these  studies,  whether  for 
analysis  of  amino  acid  sequence  or  analysis  of  post- 
translational  modification.  This  precaution  was  applied 
to  limit  sample  exposure  to  conditions  that  enhance  Asn 
deamidation,  specifically  the  high  ionic  strengths  ofgua- 

mide.®6’7®lydr0Chl°ride’  dithlothreitoI-  and  iodoacete- 

After  collection,  mass  spectral  data  were  analyzed  in 

pe°ctrTd  !°  Venfy  t0tal  amin°  aCld  Se^e—  FiJst  the 
spectral  data  were  parsed  against  a  custom  database 

containing  only  the  rPA  amino  acid  sequence.  Filters  of 
the  Spectrum  Mill  search  engine  (Agilent  Technolojes) 
ere  set  to  allow  up  to  five  missed  cleavages,  passing 
peptide  scores  greater  than  15  in  combination  with  SPI 
ored  peak  intensity)  greater  than  70%.  These  criteria 
returned  results  of  high  confidence  and  such  peptMel 
were  considered  to  be  authenticated  (“validated’^by  soft¬ 
ware  nomenclature).  In  the  second  step,  poorer  quality 
matches  of  lower  SPI  and  peptide  scores  were  analyzed 
wi  h  instrument  LC  software  (ChemStation,  Agilent) 
and  manual  inspection.  Both  precursor  and  product  ions 
were  confirmed  and  the  validity  of  the  peptide  was 
decided  based  on  manual  inspection  of  peaks  The  prod 

S  ‘hat  Sh°wed  fragment  J  J'ft 

validate!  £;SeTUOn  Were  f™*  «*  * 
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Verification  of  deamidation  and  other 
ammo  acid  modifications 


Instrument  and  system  set-up  for  collection 
of  deamidation  spectra 

An  Agilent  1100  Series  LC/MSD  Trap  SL  was 
equipped  with  an  electrospray  ionization  (ESI)  source  to 

SS  7’TS  yfb“  ^  Bo»0  »>- 

Sr  a  To"S”f'd  at  7?  “d  SSS 

B  (0  L  f  T  aCld)  and  Premixed  runnmg  buffer 

mim  28  y"1'1'!  “  acetonitriIc).  LC  method:  200  pL/ 

min,  28  pL/mjection;  gradient  from  5%  B  to  45%  B  over 

mm.  The  MSD  method:  positive  polarity  with  ICC  oa 
rameters  of  100,000  target  value,  sVoOmTi^ 
accumulafion  time  at  300-1800  m/z  scan  range  with 
v  average®’  Normal  scan  resolution  (13,000  (m/z)/ 
c),  0  psi  nebulizer;  10  psi  drying  gas,  350  C  capillary 
temperature,  Smart  Parameter  settings  of  920  m/z  tar- 

SormmrSH^ l00%  compound  stability,  140%  trap  drive 
normal  optimization;  autoMS(n)  on,  0.15-min  cycle  time’ 

hLEu7J  n- 26'793  ,hr“h0,d 

nreshold  relativity,  preoperation  SPS,  acquisition 

.  1  3°'V  fraSmentation  amplitude,  wia 

Advanced  settmgs  preference  for  doubly  charged  precur 

rirtenrsofl0grOUP  ^  °f  five’  ^t£n  Pa-' 
rameters  of  10  scan  averages  at  maximum  resolution,  5 

MS  n)  averages,  and  4.0  m/z  isolation  width  Back¬ 
ground  ions  were  excluded  from  MS(n)  analysis  as  deter¬ 
mined  empirically.  ^  eter 
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For  each  peptide  that  was  flagged  by  automated  analy¬ 
sis  to  possibly  contain  modified  amino  acids,  raw  spectral 
data  were  inspected  manually  against  an  rPA  custom 
tabase  usmg  Spectrum  Mill  entries  for  modifications 
n  an  Asp  deamidation,  succinimidyl  intermediate 

l  aZt  °n'  '7™™'  a“‘>”a«°" 

.7  ™“  1  The"”  °r  "  the  ™P»"d«»g  change 

accoTdi  to  th  fTr°Pnate  amin°  add  substitution 
dure  (T)  P!ef  ‘ng  Step-Wlse  verification  proce- 

+2  charge  f°r  precursor  Peptide  ions  of 

,  f  harge’  but  tnPle  charge  ions  were  used  when  at 
eas  wo  fragmentation  ions  were  observed  with  the 
xpected  masses  and  none  with  a  contradictor  mass-  (2) 
The  extracted  ion  chromatogram  (EIC)  for  each  of  these 
confirmed  peptides  was  plotted  by  instrument  software 
to  confirm  multiple  forms.  The  observed  precursor  ion 
using  mass  ranges  according  to  specific  modification  an" 

cerned^h6  T  w  I!™3’  Specifica%  most  analyses  con¬ 
ing  a  defta  °y  f  deamidated  P^ide  ions  hav- 
mg  a  delta  mass  of  +1  (A  m/z  =  +0.5).  The  EIC  range 

set  Si  d'“lidati''”  P™d“*  .« 

<™ ~  SSrzzs  5 

peptide  L450-R468  which  contains  one  succinimide,  dou- 
y  e  arged  peptide  ions  having  delta  mass  of  -17  (A  m/ 

?,  ~  8.5)  employed  different  EIC  limits:  range  B  - 

^precursor  mass  -16)/2  -  1  to  (precursor  mass  +1V2  ^ 

6].  Analysis  of  peptide  E704-K722  carrying  two  sites  of 
moderation  utilized  triply  charged  succinimide  peptide 
ons,  also  having  a  delta  mass  -17  (but  A  m/z  =  -5^6) 
sed  range  C  -  [(precursor  mass  -15)/3  -  1  to  (precur’ 
sor  mass  +2)/3  +  2],  The  triply  charged  ions  were  used 

intensftLasUie5^22ytfragmente<1  3nd  pr0vided  Peak 

douhTv  ,h  a  great6r  than  those  ^served  for 

P Xeslhlfd  dIOnS;  rblmg  b6tter  ^  -kulations. 

medtTe  t.  °W  evidence  of  succinimide  inter- 

ring  a  range  that  c°ntained 

ntirll  p  deamidated  peaks  and  doubly  charged 
peptides.  Every  site  was  examined  for  the  presence  of 
LC-separated  precursor  peptide  ions,  two  peaks  for 
nafive  versus  deamidated,  and  to  two  more  peaks  for  iso- 

HPLO  trT nit f  -  Differential  reversed  phase  (RP- 
HPLC)  retention  between  the  native  and  modified  states 

madSa°tionrVaeithf0r  T**  C3ITying  3  confi™ed  dea¬ 

midation,  although  not  every  deamidation  showed  iso 

ingPa°rsiSnglemErcde  *7™'  assignmenta  show- 

g  EIC  peak  were  eliminated;  and  (3)  The  v 

and  6  product  ions  of  MS/MS  spectra  for  each  precurso^ 

peptide  were  viewed  along  successive  overlapping  frag 

r;,  t7  r ,hat  thi7 Da  ahi«  2 

’  . ,  Da  mass  sbift  (succinimide),  aligned  with 
the  residue  expected  to  be  modified,  con  urrStlv  M  o 
confirming  all  surrounding  ,miiw  acid  “"enct 
quence  on  that  peptide.  Rare  isobaric  precursor  ions 
found  to  bo  peptides  of  wholly  different  a„,7no  ac“ 
sequence  were  eliminated  from  further  consideration.  Sum 
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ilar  methodology  was  applied  to  confirm  Met  oxidation 
but  using  a  delta  mass  of  +15.99  Da. 

Measurement  of  absolute  deamidation  and 
calculation  of  percent  deamidation 

Once  native  and  deamidated  residues  were  confirmed 
and  verified  in  association  with  separated  precursor 
ions,  the  absolute  amounts  of  native  and  deamidated 
species  were  estimated  by  integrating  the  peptide  EIC 
peak  with  instrument  default  settings.  Peak  areas  were 
calculated  from  the  peptide  precursor  ion  EIC  using  the 
monoisotopic  mass  with  custom  range  window  (ranges  A, 

B,  and  C  specified  above)  so  as  to  integrate  the  most 
intense  portions  of  the  isotopic  envelope. 

Relative  deamidation  at  each  confirmed  site  was  then 
calculated  from  EIC  chromatographic  peak  areas  as  a 
percentage,  comprising  the  quotient  of  the  native,  or 
deamidated  peptide  ion  in  relation  to  the  sum  of  areas 
for  all  pertinent  ions  [Eq.  (1)], 

Z)j  =  {{Ad,) /{Ad,  +A^)]100  (1) 

where  D,  is  the  percent  deamidation  calculated  at  posi¬ 
tion  i,  Ad  is  the  peak  area  of  the  precursor  ion  con¬ 
firmed  to  be  deamidated  at  position  i  by  containing  Asp 
in  place  of  the  expected  Asn,  and  AlV>  is  the  peak  area  of 
the  precursor  ion  confirmed  to  be  native  by  containing 
Asn  at  position  i,  and  all  residues  of  the  peptide  outside 
position  i  confirmed  according  to  expected  amino  acid 
sequence.  Where  observed,  AD,  also  included  EIC  peak 
areas  for  isoAsp  or  succinimide,  typically  situated  before 
the  Asn  peak  and  after  the  Asp  peak,  respectively,  due 
to  differing  retention  times  by  RP-HPLC.66  Percent  dea¬ 
midation  at  Gin  is  calculable  by  a  corresponding  equa¬ 
tion,  but  was  not  needed  as  Gin  deamidation  was  not 
observed  in  rPA.  The  same  basic  equation  was  used  to 
measure  percent  oxidation  of  Met. 

Biochemical  and  Biological  Activities 

Native  PhastGels  were  used  to  evaluate  the  ability  of 
furin  to  cleave  different  preparations  of  PAa3,  of  cut  PAe3 
to  assemble  into  (PA^L  heptamers,  and  of  the  heptamer 
to  bind  LF  in  forming  the  LFIPA^L  holoenzyme  com¬ 
plex,  as  previously  described.21  The  biological  activity  of 
various  preparations  of  rPA  to  function  in  association 
with  LF  as  anthrax  lethal  toxin,  LFlrPA^b,  was  tested 
by  measuring  the  viability  of  macrophage  cells  after  ex¬ 
posure,  in  an  improved  method  originally  reported  by 
Friedlander.74  Dilutions  of  rPA  were  prepared  in  tripli¬ 
cate  and  combined  with  40  ng/mL  of  LF  (final  concentra¬ 
tion)  in  cell  culture-treated,  96-well  Costar  plates 
(Coming  Incorporated,  Acton,  MA).  Dilutions  were  pre¬ 
pared  in  Dulbecco’s  Minimal  Essential  Medium  Com¬ 
plete,  i.e.,  with  the  following  additions:  high  glucose 
(Invitrogen-Gibco,  Gaithersburg,  MD)  containing  7.5% 
heat-inactivated  fetal  bovine  serum,  25  m M  HEPES, 
4  m M  glutamine,  100  pg  per  mL  of  streptomycin,  and 
100  U  per  mL  of  penicillin.  One-hundred  microliters  per  well 


was  then  transferred  to  cell  culture-treated  96-well 
plates  seeded  with  5  x  104  of  J774A.1  cells  in  200  pL 
volumes  18-22  h  before  testing  after  removing  the 
culture  medium.  Cells  were  maintained  in  Dulbecco  s 
Minimal  Essential  Medium  Complete  without  25  m M 
HEPES.  After  incubating  the  plates  for  4  h  in  a  humidi¬ 
fied  incubator  set  at  37°C  and  5%  C02,  25  pL  per  well  of 

3-[4,5-dimethylthiazol-2-yl]2,5-diphenyltetrazolium  bro¬ 
mide  (MTT)  at  5  mg/mL  of  PBS  was  added  and  the 
plates  reincubated  for  2  h.  Cells  were  then  lysed  and  the 
insoluble  reduced  formazan  was  made  soluble  by  adding 
100  pL  per  well  of  lysing-solubilization  buffer.75  After  an 
overnight  incubation,  the  absorbance  values  at  570-690 
nm  for  each  well  were  collected  by  using  a  BioTek  312e 
microplate  reader  (BioTek  Instruments,  Winooski,  VT). 
The  protein  concentration  of  each  PA  test  sample  that 
resulted  in  50%  killing  of  cells  in  the  presence  of  LF  was 
calculated  (EC50)  using  the  4-parameter  logistic  regres¬ 
sion  analysis  of  SoftMax  Pro  (Molecular  Devices,  Sunny¬ 
vale,  CA). 

Computational  Methods 

Predicted  rates  of  Asn  deamidation  in  the  monomer 
were  calculated  using  the  algorithm  of  Robison  and  Rob¬ 
inson  available  at  www.deanndation.org.  Correlations 
between  percent  deamidation  and  negative  charge  or 
biological  activity  were  calculated  using  the  Pearson 
product-moment  correlation  in  SAS  software,  Version  9.1 
(SAS  Institute,  Cary,  NC).  Modified  amino  acids  of  PA 
monomer  and  heptameric  PA  prepore  (PDB  code  1ACC 
and  1TZO)  were  viewed  and  studied  using  Open-Source 
PyMOL™  Version  0.97  (www.pymol.org)  (DeLano  Scien¬ 
tific  LLC,  SanCarlos,  CA). 

RESULTS  AND  DISCUSSION 
Visual  Analysis  of  PA  Isoforms  Indicates  Rapid 
Protein  Modification  in  vivo 

Almost  two  decades  ago,  Leppla4  reported  the  exis¬ 
tence  of  natural  PA  protein  variants  during  the  first 
isolation  and  characterization  of  B.  anthracis  exotoxin 
proteins.  PA  heterogeneity  appeared  as  micro-charge  iso¬ 
forms  of  incrementally  greater  negative  charge,  as  first 
detected  by  ion  exchange  chromatography  and  nondena¬ 
turing  slab  gel  electrophoresis.76  These  isoforms  were 
subsequently  characterized  in  detail  with  CEIF  (Fig.  2), 
native  PhastGel23  (Fig.  3),  and  2-DE  (Fig.  4).  Because 
isoforms  of  natural  and  rPA  were  all  seen  to  share  the 
expected  Mr  of  83  kDa,4’23'70  it  was  suggested  that  heter¬ 
ogeneity  could  be  caused  by  post-translational  modifica¬ 
tion,  most  fittingly  deamidation.23  Although  supportive 
evidence  for  deamidation  was  recently  reported,23’  this 
hypothesis  was  difficult  to  prove  due  to  the  small  associ¬ 
ated  delta  mass  and  was  not  definitively  and  specifically 
demonstrated  for  PA  protein  until  this  work,  as 
described  herein.  First,  we  considered  how,  in  light  of 
certain  evidence  for  deamidation,  the  evaluation  of  visi¬ 
ble  isoforms  revealed  two  unexpected  details  of  rPA 
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pi9  ■  2.  Fractionation  ot  rPA  isoform*;  hw  pirr  ,  . 

R  &  D  lot  of  purified  rPA  nrntein  Z*  CIER  lso,orm  classes  in  an 
pi  values,  calculated  as  described  in  c*  6  3S  dis,inct  Peaks-  Measured 
three  largest  peaks  are  markedbj  arrowrS'sland^^68'  I?'  ,he 
y  intercept  =  10.24  sIodp  -  ~n  nnoo  „fet'  standard  curve  having 
and  pis  of  analytes  used  in  rnnctr  *•  3r"l^  ~  °  "28-  The  mobilities 

drase  isoform  dI  5  9  c  mi_.  „  ,  1 1  P  y-5’  1 1  m,n.  carbonic  anhy- 
24.3  min;  CCK  inking  pep, id™"  pil^™  em^8^0'0™'  P'  54 

5“  sample  ^Sy^ 

of  purged  rPA  contained  substantial  modification 
Fractionation  by  CEIF  of  purified  R  &  D  grade  rPA 
protein  showed  four  clear  chromatographic  -aka  Zo 
dominant  and  two  minor  (see  Fig.  2).  However  peak 
skewing  and  shouldering  indicated  incomplete  separa 
ion  and  the  probable  existence  of  additional  minor  coe- 

S"Tf  Unr6S0l-d  sPecie®-  The  three  most  abundant 
CEIF  isoforms  each  differed  in  pi  bv  ~0  1  nH  ,, 

all  values  were  Cose  to  (within  abS^pH  un£  but  meTs 

T^r  r  y  ptSS’  than  the  calculated  Pi  of  5.88  for  rPA  71 
Nat!ve  PhastGel  analysis  was  developed  as  a  convenient 
and  robust  assay  of  protein  purity  and  aualitv  in 
production  of  pharmaceutical  grade  rPA  protein  manu¬ 
factured  under  GMP  controls  for  testing  rPA  as  tlTnext 
anthrax  vaccme  antigen.232470  Method  optimization 
through  process  development  finally  reduced  the  number 
and  apparent  mtensitie8  of  ViSlble  charge-class  species 
(isoform  complexity)  in  rPA  to  two  principal  isoforms 
with  insubstantial  contamination.  The  single  most  effec 
five  procedure  used  to  reduce  isoform  compkSy  was 
the  introduction  of  an  early  ammonium  sulfate  cut  as  a 

This^iirtafied  th*  ^  Unpubllshed  results). 

for!^«  w  d  he  aCCUmulation  of  additional  acidic 
forms  visible  m  various  R  &  D  lots,  and  effectively  con- 

ntrated  the  two  isoforms  which  were  seen  to  be  univer- 
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E)  and  aged  (B,  D;;> LC-MS/MS  for  deamidation  at  N537  as  contained  on  peptide  529-1 AFGFNEPnGNL- 
QYQGK9545  are  indicted  by  histogram  bars:  dark,  deamidated;  light,  native.  Relative  ion  peak  areas  (ordinate)  are  plotted  as  a  function  of  the  iso- 
form  spot  analyzed  (abscissa). 


animal  vaccination  and  challenge  studies.23  All  of  these 
results  clearly  demonstrated  that  isoform  complexity 
increased  with  protein  degradation,  decreased  with  pro¬ 
tein  purity,  and  began  in  vivo  early  during  expression. 
Importantly,  each  of  these  techniques  revealed  the  pres¬ 
ence  of  minor  rPA  species  of  less  negative  charge  than 
the  components  of  highest  abundance,  evident  as  small 
peaks  preceding  the  major  CIEF  peaks  (Fig.  2),  minor 
bands  above  the  major  PhastGel  bands  (see  Fig.  3),  and 
minor  spots  of  higher  positive  charge  than  the  major  2- 
DE  spots  [Fig.  4,  e.g.,  panel  (D)  spots  7  and  8].  Review  of 
rPA  molecular  charge  heterogeneity  by  these  three  frac¬ 
tionation  techniques,  and  by  preparative  ion  exchange 
chromatography  reported  earlier,23,70  confirmed  that  the 
primary  forms  of  purified  PA  were  not  the  isoform  spe¬ 
cies  of  least  net  negative  charge,  and  consequently,  were 
not  the  least  modified  form  of  PA.  Thus,  by  the  following 
reasoning,  we  conclude  that  most  GMP  grade  PA  protein 


is  modified  in  vivo,  and  largely  to  the  extent  apparent 
in  vitro  after  purification. 

First,  the  two  principal  isoforms  predominated  in  all 
preparations  but  did  not  comprise  the  entire  content  of 
any  purified  rPA  preparation.  These  isoforms  were  seen 
to  dominate  rPA  isoform  populations  from  the  earliest 
point  of  detection  during  fermentation,23’24'70  and  for  this 
reason  they  were  previously  designated  “isoform  1”  and 
“isoform  2.”23  While  these  two  isoforms  together  comprise 
99.72%  of  the  total  protein  in  GMP  grade  rPA,24  and 
likely  as  much  of  high  purity  rPA  produced  by  a  commer¬ 
cial  vendor  [List  Biologicals,  Fig.  6(A)],  they  were  also 
accompanied  by  very  faint,  less  acidic  isoforms  as  was 
apparent  in  freshly  thawed  GMP  grade  protein  (Figs.  3 
and  4),  aged  GMP  protein  [Fig.  5(A)],  rPA  produced  by  a 
commercial  vendor  [Fig.  6(A  sample  a)],  as  well  as  R  &  D 
lots  of  rPA  [Figs.  2,  6(A  samples  b-d].  Because  isoforms 
ostensibly  represent  incremental  post-translational  modi- 
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Fig  5.  Multiple-site  deamidation  and  decreased  biological  activit 
correlate  with  increased  isoform  complexity  of  aged  GMP  rPA.  Compar 
son  of  isoform  complexity  and  biological  activity  with  percent  deamida 
tion  at  seven  Asn  sites  in  freshly  thawed  versus  aged  GMP  rPA  (A 

tpfn'VRPhfStGelS,  °f  Wh°le  GMP  ^  Protein;  Left.  freshly  thawed  prc 
tein,  Right,  accelerated  aged  protein.  Horizontal  guide  lines  aide  ii 
comparing  vertical  positions  of  protein  isoforms  between  gels.  (B)  Bio 
togical  activity  measured  as  ECso  (ng/mL)  with  standard  deviation  ii 
brackets.  (C)  Relative  deamidation  measured  in  each  protein  sampl. 
(abscssa)  expressed  as  percent  deamidated  (ordinate).  Samples,  verb 
Mby  aligned  with  gel  images  in  A,  comprise:  a,  freshly  thawed  GMF 

M^b’Ae5ed  °r  acce,eratecl  aging.  Residues:  1,  N602;  2  N408'  3 
N719;  4,  N601;  5,  N466;  6,  N713;  7,  N537;  M,  mean. 


fixations  that  impart  progressively  greater  negative 
charge,  the  less  negative  isoforms  were  presumably  less 
modified  (i.e.,  more  native)  than  the  dominant  and  highly 
visible  isoforms.  Second,  the  uppermost  prominent  iso¬ 
form  m  all  final  preparations  coaligned  by  native  Phast- 
Gel,  as  seen  by  comparing  isoform  distribution  among  R 
&  D  lots,  the  GMP  rPA  lot,  and  a  lot  produced  by  a  pro¬ 
prietary  and  possibly  different  process  [Figs.  5(A)  and 
6(A)],  Therefore,  these  prominent  bands  had  the  same 
net  molecular  charge.  Third,  this  same  predominant 
upper  isoform  band  [band  3  in  Fig.  3(A)]  aligned  with  the 
upper  most  prominent  band  visible  in  the  earliest  stages 
of  expression  sampled  during  fermentation,  previously 
called  “isoform  1”  and  shown  in  Figure  1  of  Ribot  et  al.23 
Therefore,  the  putative  modification  causing  this  isoform 
occurred  naturally  and  rapidly  during  expression  in  the 
bacterium.  Fourth,  this  same  uppermost  dominant  iso¬ 
form  of  rPA  had  a  measured  pi  value  notably  less  than 
that  calculated  in  silico  (5.49  vs.  5.88),  was  preceded  by 
two  minor  isoforms  of  observed  pi  closer  to  the  calculated 


rPA  sample  (lot) 

,6  ^ultiPle-site  deamidation  and  decreased  biological  activity 
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N602;  2,  N408;  3,  N719;  4,  N601;  5,  N466;  6,  N713;  7,  N537;  M,  mean.  ’ 


value,  and  was  succeeded  by  isoforms  with  increasingly 
acidic  pi  values  (Figs.  2  and  3).  Thus,  most  of  the  PA  pro¬ 
tein  produced  during  expression  and  present  after  purifi¬ 
cation  exhibited  charge  states  that  were  more  acidic  (neg¬ 
atively  charged)  than  those  calculated  from  amino  acid 
sequence  and  were  therefore  more  modified  than  nascent 
rPA  protein.  Finally,  in  every  preparation  we  studied,  the 
trend  of  increasing  isoform  complexity  comprised  an 
acidic  shift  in  isoform  distribution  beyond  that  of  the 
principal  isoforms,  whether  observed  during  bacterial 
p-owth,  collected  as  fractions  across  a  primary  peak  dur¬ 
ing  isoform  purification,23  created  inadvertently  through 
andlmg  or  extended  4°C  hold  during  process  develop- 
ment  [Fig.  6(A)],  or  intentionally  produced  by  freeze/ 
thaw/hold  cycling  to  achieve  accelerated  aging  [Fig.  5(A)], 
We  conclude  that  PA  protein  was  rapidly  and  naturally 
converted  to  several  successively  acidic  states  in  vivo 
two  of  which  predominated  and  were  enriched  during  pu¬ 
rification.  All  of  these  findings  are  concordant  with  the 
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earlier  suggestion  that  the  probable  cause  of  PA  hetero¬ 
geneity  is  multiple-site,  amino  acid  deamidation.23  Two 
opposing  explanations  concur  with  these  findings:  either 
the  principal  components  were  not  fully  native  and  nega¬ 
tively  charged  modifications  were  the  only  alteration 
present,  as  proposed  here,  or  there  may  have  also  been 
hypothetical  modifications  imparting  positive  charge.  We 
offer  direct  and  definitive  proof  for  the  first  hypothesis. 
Such  evidence  required  the  development  of  a  robust  LC- 
MS/MS  method  for  detecting  and  inspecting  every  site  of 
deamidation  with  an  empirical  approach  for  their  quali¬ 
tative  and  semiquantitative  analysis.  With  these  tools, 
gel-isolated  isoforms  were  inspected  for  variation  in  dea¬ 
midation  in  relation  to  isoform  complexity  to  define  the 
chemical  source  of  PA  heterogeneity. 

Complete  Amino  Acid  Sequencing  of  rPA  Protein 
by  LC-MS/MS 

Evaluation  for  deamidation  began  with  the  complete 
sequencing  of  rPA  by  MS.  Analysis  of  MS/MS  spectral 
data  for  GMP  protein  digested  with  trypsin,  chymotryp- 
sin,  Glu-C  (V8  protease),  or  combinations  of  these  three 
proteases,  and  parsed  automatically  against  a  custom 
rPA  database  confirmed  86%  of  the  predicted  amino  acid 
sequence.  Four  expected  peptides  were  not  observed 
within  set  limits  of  the  mass  analysis  program  including 
22-YSFDLNFQAPMWTSSTTGDLSIPSSELENIPSENQY- 
FQSATWSGFI-KVK,  111-GRLYQIKIQY,  133-K,  and  158- 
QKSSNSRKK.  However,  we  were  able  to  detect  and  frag¬ 
ment  these  peptides  and  extend  amino  acid  sequence 
coverage  to  100%  by  eliminating  the  peptide  ion  match¬ 
ing  tolerance  limits,  instead  using  instrument  software 
and  manual  inspection  to  confirm  these  ions  (see  Materials 
and  Methods).  Merging  of  all  data  resulted  in  total 
amino  acid  sequence  confirmation  with  the  observation 
of  every  possible  source  for  deamidation,  as  shown  in 
Figure  7.  All  68  Asn  residues  and  all  30  Gin  residues 
were  observed,  and  those  showing  mass  deviation  were 
inspected  for  confirmation  as  being  either  native  or  dea¬ 
midated  through  redundant  mass-based  amino  acid 
sequence  analysis  with  manual  verification.  Automated 
analysis  of  the  spectral  data  flagged  61  Asn  and  six  Gin 
for  potential  deamidation,  plus  one  Met  for  potential  oxi¬ 
dation  (Fig.  7).  However,  most  of  these  flagged  modifica¬ 
tions  were  erroneous  assignments,  as  determined  by  the 
verification  process  discussed  next.  No  deletions  were 
detected,  and  a  directed  search  for  other  possible  amine- 
specific  modifications  listed  in  the  search  engines 
revealed  no  other  post-translational  modifications. 

A  Simple  Method  to  Verify  Deamidation: 
Demonstration  of  Site  Specific  and  Multiple 
Asn  Deamidation  in  PA 

While  the  previously  reported  enzymatic  detection  of 
isoAsp  provided  strong  supportive  evidence  of  deamida¬ 
tion  in  PA,23  it  did  not  confirm  the  existence  of  deamida¬ 
tion  or  reveal  its  nature.  Alternatively,  MS  techniques 
such  as  MALDI-TOF-MS  that  do  not  preseparate  peptide 


EVKQERRLLME3E3SSQGLl.GTTE3DLHEQEPHWTS3TTG0L3IPS3El.E*IPSERQTEQOAIW3GriKViat3BETTr4  BO 
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■  HOW.AEUm»irrVLPEIKIJEKM3LIRPERrgTPRiw>IAVGADESWKEAHREVI»33TEGLI.miPKDIRKILS  840 

"  _  713  718 
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Fig.  7.  Summary  of  LC-MS/MS  analysis  of  rPA  for  amino  acid 
sequence  and  post-translational  modifications.  Shown  is  the  amino  acid 
sequence  of  rPA  (accession  PI  3423)  as  experimentally  verified  by  LC- 
MS/MS  in  universal  single  letter  code,  beginning  with  the  N-terminal 
amino  acid  remaining  after  natural  removal  of  the  leader  sequence  (i.e., 
PAaa).  Special  nomenclature  and  markings  as  follows:  n,  Asn  residue, 
identified  by  coordinate,  experimentally  confirmed  to  be  partially  deami¬ 
dated  by  also  containing  Asp;  m,  residue  observed  to  be  partially  oxi¬ 
dized  by  also  containing  oxyMet;  bold  type,  amino  acid  sites  predicted 

to  be  modified  by  automated  MS/MS  analysis;  solid  underscore  ( _ ), 

tryptic  peptides  containing  7  confirmed  sites  of  deamidation;  superscript 
number  and  asterisk  (*),  residue  coordinate  confirmed  to  be  modified 
as  described  in  Experimental  Procedures.  Findings  are  compiled  from 
data  collected  and  analyzed  on  fragments  from  digestion  with  trypsin, 
chymotrypsin,  or  endoproteinase  Glu-C  (V8  protease)  of  GMP  rPA  as 
described. 

forms  by  LC  or  analyze  ion  composition  by  MS/MS  frag¬ 
mentation  can  confirm  deamidation  but  cannot  discrimi¬ 
nate  the  location,  nature,  or  quality  of  the  putative 
amino  acid  deamidation.  Calculation  of  deamidation  by 
such  methods  rests  on  comparing  isotope  peak  abundan¬ 
ces  between  the  monoisotopic  A  peak  and  all  major 
peaks  of  the  isotopic  series,  which  principally  includes 
the  A+l,  A+2,  A+3  peaks  for  average-sized  peptides. 
However,  this  measurement  of  peak  abundances  is  only 
provisional  without  better  than  5  millimass  unit  accu¬ 
racy  or  certain  knowledge  of  elemental  composition.  This 
is  because  the  A+l  ion  of  the  12C  form  for  a  truly  deami¬ 
dated  peptide  (A  m  =  +0.9841  Da  for  a  single  Asn  dea¬ 
midation)  and  the  first  heavy  isotopic  form  of  the  native 
peptide  (e.g.,  A  m  =  +1.003  Da  for  a  single  13C)  are 
collected  together  within  the  same  m/z  range  window 
during  a  time  of  flight  mass  spectrometry  experiment. 
Without  high  mass  accuracy  measurement,  the  resulting 
over-representation  in  A+l  peak  area  must  be  assumed 
as  constant  between  differentially  treated  samples. 
Moreover,  this  method  cannot  distinguish  multiple  dea¬ 
midation  present  on  the  same  peptide. 

In  contrast,  we  optimized  a  simple  LC-MS/MS  method 
to  definitively  ascribe  an  observed  +1  Da  mass  shift  to 
deamidation  at  a  specific  amino  acid  residue,  even  using 
an  instrument  of  0.2  m/z  nominal  mass  accuracy.  This 
method  of  verification  was  based  on  satisfying  three  cri¬ 
teria  for  every  deamidation  while  inspecting  raw  spectra: 
criterion  (1)  the  presence  of  resolved  RP-HPLC  peaks  for 
all  variants  of  the  deamidated  peptide  in  question,  with 
delta  mass  of  +1  Da  (Asp,  isoAsp,  Glu)  or  -17  Da  (succi- 
nimide);  criterion  (2)  confirmation  among  fragmentation 
ions  that  each  mass  shift  aligned  with  the  suspected 
modified  residue  in  the  observed  peptide  amino  acid 
sequence;  and  criterion  (3)  confirmation  that  the  monoi- 
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A+l  peak  in  assigning  its  monoisotopic  mass.  Inspection 
of  the  single  flagged  Met  (M521)  using  a  mass  shift  of 
+15  99  Da  for  oxy-Met  sulfoxide  confirmed  it  was  par¬ 
tially  oxidized,  with  1.43%  oxyMet  content  measured  in 
fresh  protein  (Table  I).  The  amount  of  oxyMet  did  not 

increase  with  treatment  (Table  I). 

In  summary,  we  found  that  GMP  grade  rPA  protem  con¬ 
tains  seven  Asn  residues  with  detectable  and  quantifiable 
deamidation  (Fig.  7),  plus  an  eighth  site  just  barely  detect¬ 
able  but  too  low  to  quantify.  The  only  additionally  observed 
modification  was  a  slight  oxidation  at  a  single  Met  and 
oxyMet  does  not  impart  a  charge  to  the  protein.  These 
data  provided  definitive,  comprehensive,  and  specific 
chemical  support  for  the  long  speculated  but  unproven  hy¬ 
pothesis  that  PA  protein  isoforms  are  caused  by  ammo 
acid  deamidation.  Interestingly,  none  of  these  seven  sites 
was  seen  to  be  fully  modified,  and  each  varied  m  content 
according  to  isoform  complexity,  as  discussed  below. 


Demonstration  of  absolute  and  relative 
deamidation:  cataloging  fractional 
deamidation  in  purified  rPA  protein 

After  establishing  all  sites  of  modification  and  observ¬ 
ing  only  partial  levels  of  deamidation  in  every  case,  we 
explored  whether  these  data  could  provide  an  empirical 
measure  of  absolute  deamidation.  While  inspecting  raw 
LC-MS  and  MS/MS  data  to  verify  site  specific  deamida¬ 
tion  we  observed  that  most  peptide  variants  were  suffi¬ 
ciently  separated  by  LC  before  mass  analysis,  and  so  it 
became  evident  that  the  amount  of  deamidation  could  be 
estimated  by  direct  integration  of  EIC  chromatographic 
peak  areas  instead  of  the  typical  method  of  calculating 
differences  between  observed  and  inferred  isotope  pea 
profiles.  The  absolute  amount  of  each  peptide  form  was 
measured  independent  of  its  comparative  variant  of  in¬ 
terest,  and  relative  levels  of  deamidation  were  calculated 
directly  from  these  empirical  values.  Thus,  this  method 
circumvented  common  calculations  that  gauge  compara¬ 
tive  peak  areas  derived  from  theoretical  isotopic  profiles 
or  from  differentially  treated  samples.61'63  Furthermore, 
because  RP-HPLC  separates  peptides  according  to 
hydrophobic  character  and  not  mass  [criterion  1,  illus¬ 
trated  in  Fig.  8(A)],  observed  peptide  isoforms  of  +1 
delta  mass  comprised  only  the  deamidated  species  and 
not  a  heavy  version  (e.g„  13C)  of  the  native  species  [crite¬ 
rion  3,  illustrated  in  Fig.  8D-F],  In  this  method,  isotope 
analysis  was  used  only  to  support  the  interpretation  of 
MS/MS  fragmentation  data  in  confirming  the  identity, 
purity,  and  correct  mass  assignment  of  each  peptide  ion. 
Calculations  across  all  confirmed  sites  of  modification  m 
rPA  confirmed  differing  levels  of  partial  deamidation, 
and  in  no  case  was  a  native  peak  absent  (i.e.,  100%  dea¬ 
midation),  including  peptide  529-545,  which  carried  the 
Asn  of  highest  observed  deamidation  [Figs.  5(C)  and 
6(C)].  The  clarity  of  measuring  deamidation  in  this  man¬ 
ner  is  visually  illustrated  in  Figure  8  for  this  same  pep¬ 
tide,  wherein  freshly  thawed  GMP  grade  protein  had  a 
larger  native  N537  peak  (peptide  isoform  containing  only 


Asn  at  position  537)  than  did  treated  protein,  while  the 
treated  protein  had  a  larger  deamidated  peak  (peptide 
isoform  containing  only  Asp  at  position  537)  than  di 
fresh  protein.  As  calculated  from  only  confirmed  peak 
areas  using  Eq.  (1),  the  measured  percent  deamidation 
at  this  site  was  35.5%  in  fresh  protein  and  85.2%  in 
treated  protein  [Fig.  5(0,  Table  I].  The  presence  of  both 
native  and  deamidated  forms  of  a  peptide  separated  by 
charge  via  LC  before  mass  analysis  [Fig.  8(A)]  was 
observed  for  every  proteolytic  fragment  carrying  a  con¬ 
firmed  Asn  deamidation.  However,  before  measuring 
each  site  for  absolute  and  relative  deamidation,  we  first 
verified  that  these  calculations  were  reproducible  using 
model  synthetic  peptides  of  differentially  defined  compo¬ 
sition  as  follows. 

While  electrospray  ionization  MS  can  provide  quanti¬ 
tative  measures  of  ions  within  certain  density  ranges  o 
the  instrument  ion  trap,  we  nevertheless  tested  the 
hypothesis  that  the  proposed  method  could  accurately 
estimate  levels  of  deamidation.  For  this  purpose,  we  per¬ 
formed  measurements  on  synthetic  peptides  that  mim¬ 
icked  the  tryptic  fragment  containing  the  site  of  highes 
observed  deamidation  (529-IAFGFNEP1VGNLQYQGK- 
545),  such  that  one  peptide  contained  the  native  form 
(N537)  while  the  second  contained  the  fully  deamidated 
and  PIMT-corrected  form  (D537).  Dilutions  and  measure¬ 
ments  of  each  peptide  in  isolation  confirmed  that  a 
standard  curve  could  be  constructed  for  each  peptide 
with  good  fit  to  linearity  (native  N537  peptide:  R  - 
0.9975;  deamidated  D537  peptide  R2  =  0.9918).  This  lin¬ 
ear  range  encompassed  most  of  the  ion  intensities 
observed  experimentally  for  rPA-digested  peptides  [Figs. 

3  4,  5(C),  6(C),  Table  I].  Moreover,  mixing  the  two  syn¬ 
thetic  model  peptides  in  different  ratios  and  measuring 
them  as  “unknown  samples”  demonstrated  that  their 
measured  fractional  compositions  aligned  with  calculated 
stoichiometric  compositions  (data  not  shown).  These 
results  confirm  the  proposal  that,  at  least  for  these  two 
peptide  mimetics  of  native  and  deamidated  tryptic  frag¬ 
ment  529-545,  this  method  dependably  and  accurately 
measured  the  absolute  and  fractional  amounts  of  deami¬ 
dation  at  residue  N537.  Without  specific  demonstration, 
we  propose  that  correspondingly  similar  representation 
can  be  extended  to  other  rPA  peptides  carrying  sites  of 
deamidation.  Thus,  through  direct  measurements  of  sep¬ 
arated  native  and  deamidated  peptide  peak  areas,  this 
method  provided  an  estimate  of  absolute  and  relative 
deamidation  at  every  observed  modified  position  without 
dependence  on  deduced  or  differential  isotope  profiles. 

Absolute  and  relative  levels  of  deamidation  were  meas¬ 
ured  across  all  modified  Asn  residues  in  freshly  thawed 
GMP  grade  rPA  protein.  All  Asn  deamidations  reacted  to 
a  different  extent  and  could  be  arranged  in  the  follovnng 
rank  order  of  net  conversion  to  Asp:  N537  >  N713  > 
N466  >  N719  >  N601  >  N408  >  N602  [Table  I,  Fig. 
5(C)].  Using  this  single  point  of  observed  Asp  content  as 
an  approximate  measure  of  an  Asn  residue’s  propensity 
to  undergo  deamidation,  these  data  appear  to  differ  from 
the  rank  order  for  rates  of  deamidation  predicted  by  the 
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Robinson  algorithm  [Fig.  10(A)],  In  this  regard,  Zomber 
et  al.63  recently  reported  tantalizing  MALDI-TOF-MS 
fingerprinting  data  for  PA  that  concur  with  and  differ 
from  our  findings  of  Asp  deamidation,  first  reported  ear¬ 
lier  and  fully  described  here.  In  agreement,  deamidation 
at  N713  and  N719  were  inferred  as  the  cause  of  a  +2  Da 
shift  of  peptide  704-723  which  contains  four  Asn  resi¬ 
dues  that  could  possibly  degrade.  However,  without  MS/ 
MS  sequencing,  individual  levels  of  deamidation  and 
direct  proof  of  which  Asn  had  actually  modified  could 
not  be  ascertained  directly  on  digested  protein,  but 
instead  was  deduced  from  theoretical  predictions  of  dea¬ 
midation  and  from  measures  of  model  synthetic  pepti¬ 
des.38  Aged  synthetic  peptides  were  reported  to  exhibit 
100%  deamidation  at  these  two  sites,63  where  we  found 
no  more  than  37%  deamidation  at  either  of  these  sites  in 
treated  PA  protein  of  high  isoform  content  (Table  I,  Figs. 

5  and  6).  Also  in  concordance  with  data  reported  here, 
no  evidence  was  found  for  deamidation  at  N10  or  N306, 
which  is  in  sharp  contrast  to  predictions  from  the  Robin¬ 
son  algorithm  that  calculates  both  of  these  sites  to  be 
among  the  most  reactive  Asn  residues  in  PA  (Fig.  10). 
Our  findings  are  also  comparable  with  respect  to  possi¬ 
ble  deamidation  at  N162,  a  site  of  possible  significance 
as  it  juxtaposes  the  furin  cleavage  site.  Namely,  neither 
approach  found  unequivocal  evidence  of  deamidation  at 
this  position  within  fresh  PA  protein.  Using  aged  isolates 
of  cleaved  and  isolated  PA20,  an  observed  +1  Da  mass 
shift  for  peptide  149-167  was  attributed  entirely  to  one 
(N162)  of  two  Asn  residues  present  on  the  peptide,  based 
on  theoretical  rates  of  deamidation.63  Instead,  inhibition 
of  furin  activity  was  demonstrated  by  substitutions  anal¬ 
ogous  to  position  N162  of  synthetic  peptide  mimetics. 

We  found  no  convincing  evidence  of  modification  at  N162 
even  in  aged  PAgs  protein  (Fig.  7),  and  therefore  propose 
that  the  observed  inhibition  of  furin  activity  for 
degraded  PA  protein  may  be  due  to  significant  deamida¬ 
tion  at  other  positions,  as  described  below.  Thus,  find¬ 
ings  reported  previously  and  herein  reflect  two  areas  of 
agreement.  First,  that  the  specific  locations  of  multiple 
deamidation  and  their  individual  contributions  to  an 
observed  mass  shift  in  a  protein  cannot  be  directly  and 
unequivocally  determined  for  proteolyzed  protein  using 
MALDI-TOF-MS  fingerprinting  profiles  alone.  We  show 
instead  that  such  determination  could  be  directly 
achieved  by  LC  (i.e.,  nonmass-based)  peptide  separation 
before  MS  analysis  with  concurrent  MS/MS  sequence 
analysis.  Second,  that  the  Robinson  algorithm  for  pre¬ 
dicting  deamidation  does  not  fit  experimental  data  for  at 
least  two  Asn  residues  in  PA  protein  (N10,  N306),  and 
our  own  findings  indicate  that  deamidation  at  several 
other  Asn  residues  also  do  not  conform  to  predictions,  as 
discussed  further  below. 

We  note  two  aspects  of  this  method  that  are  not  fully 
optimized  and  may  generate  inaccuracy  in  measurement 
of  deamidation  depending  on  the  particular  peptide 
being  inspected.  First,  over  estimation  of  natural  deami¬ 
dation  may  arise  from  conditions  of  treatment  in  vit™ 
that  can  enhance  deamidation  including  proteolysis. 


However,  unlike  other  published  work,  these  proteoly¬ 
sis  procedures  did  not  expose  samples  to  denaturation, 
cysteine  reduction,  or  alkylation.  The  necessary  reac¬ 
tants  (chaotrope,  reducing  agent,  alkylating  agent)  have 
high  dielectric  constants  which  is  one  factor  known  to 
enhance  Asp  deamidation  in  vitro.  Chemical  environ¬ 
ments  of  high  pH,  high  ionic  strength,  or  high  dielectric 
field  strength  stabilize  the  protonated  state  of  the 
attacking  peptide  bond  nitrogen,  which  is  the  nucleo¬ 
philic  anion  that  initiates  formation  of  the  succinimide 
intermediate  and  rate-limiting  step  of  Asp  deamidation 
(see  Fig.  1). 77,78  Variable  increases  in  deamidation  were 
observed  for  test  samples  comparing  treatments  of  4  h 
versus  the  standard  18  h  (data  not  shown).  Nonetheless, 
induced  amounts  would  cancel  one  another  upon  calcu¬ 
lating  changes  in  deamidation  since  all  proteolysis  reac¬ 
tions  used  identical  solvent,  temperature,  and  duration. 
Secondly,  under  representation  of  natural  deamidation 
may  arise  from  unaccounted  forms  of  deamidation,  spe¬ 
cifically  the  isoAsp  end  product  and  the  succinimide 
intermediate.  We  observed  these  forms  only  in  three  res¬ 
idues  of  rPA  (N466,  N713,  and  N719).  Yet  formation  of 
these  components  is  variable  and  unpredictable,  and 
their  observation  by  RP-HPLC  can  be  strongly  depend¬ 
ent  on  amino  acid  sequence  which  differs  between 
distinct  peptides.  Indeed,  their  irregular  appearance  is 
neither  unexpected  nor  unprecedented,66  and  similar 
errors  afflict  many  MS  measurements  of  protein  deami¬ 
dation.  Alternatively,  direct  and  definitive  assignment  of 
deamidation  and  isoAsp  content  can  be  achieved  using 
instruments  of  high-mass  accuracy.60,62,79  Supportive 
evidence  for  deamidation  and  isoaspartate  at  these  sites 
has  been  observed  by  electron  capture  dissociation  on  an 
FTMS  instrument  (P.  O’Connor,  unpublished  data),  as 
will  be  described  separately. 

Analysis  of  Gel-Isolated  rPA  Isoforms:  Percent 
Deamidation  Associates  With  Isoform  Net 
Negative  Charge  and  Isoform  Complexity 

Next,  we  studied  purified  isoforms  to  look  for  evidence 
that  might  support  the  supposed  correlation  between 
site  specific  deamidation  and  an  individual  isoform  or 
pattern  of  isoforms.  This  cannot  be  evaluated  using  puri¬ 
fied  rPA  protein  because  such  samples  are  demonstrably 
heterogenous,  and  a  measurement  of  deamidation  in  fact 
represents  a  population  average  of  all  individual  iso- 
forms  present  in  the  analyzed  sample.  Therefore,  iso¬ 
lated  protein  isoforms  were  excised  from  native  PhastGel 
bands  and  from  2-DE  gel  spots  then  analyzed  for  their 
individual  levels  of  deamidation.  As  in  prior  findings 
using  averaged  measurement  of  unfractionated  protein, 
Asn  N537  showed  the  highest  level  of  measured  deami¬ 
dation  in  nearly  every  isolated  charge  isoform  (Figs.  3 
and  4).  Although  the  absolute  amount  of  peptide  529- 
545  varied  due  to  differing  protein  concentration  in  each 
isoform  band,  the  percent  deamidation  at  N537  clearly 
increased  with  increasing  molecular  net  negative  charge 
of  isoform  species.  This  trend  was  evident  among  isoform 
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two  sites  of  deamidation  (N713,  N719)  and  theoretically 
could  have  contained  double  deamidation  on  a  single 
peptide.  We  cannot  exclude  the  possibility  that  such  a 
peptide  may  be  present  at  very  low  levels.  Although  four 
isoform  peptide  peaks  for  such  precursor  ions  were 
observed,  they  overlapped  substantially  and  their  El 
peak  areas  were  very  small  and  nearly  immeasurable 
using  doubly  charged  ions.  In  contrast,  tnply  charged 
precursor  ions  yielded  substantially  more  intense  and 
measurable  peak  areas,  and  were  instead  used  for  calcu¬ 
lating  deamidation  levels  (Materials  and  Methods,  Table 
I).  Moreover,  the  detection,  confirmation,  and  use  ol  suc- 
cinimide  intermediates  assisted  measurements  since  the 
delta  mass  of  -17  circumvented  confusion  otherwise 
caused  by  overlapping  heavy  isotopes  and  the  delta  mass 
of  +1  for  deamidation.  In  contrast,  peptide  601-613  also 
contained  two  deamidated  Asn  (N601,  N602)  but  three 
forms  were  well  separated  by  RP-HPLC.  Here,  differing 
amino  acid  sequences  likely  had  a  large  impact  on  the 
fractionation  and  discrimination  of  smaller  peptides.  The 
absence  of  simultaneous,  tandem  Asn  deamidation  at 
N601  and  N602  concurs  with  the  Robinson  algorithm 
which  calculates  very  low  deamidation  for  any  such  an 
arrangement.60  As  discussed  below,  most  of  the  rates 
inferred  from  measured  deamidation  in  rPA  did  not  con¬ 
cur  with  predicted  relative  order. 

In  summary,  Asn  deamidation  generally  increased  in 
association  with  rPA  isoform  complexity,  and  the  method 
reported  here  will  enable  global  inspection  of  site-specific 
deamidation  in  PA  as  a  function  of  many  treatments^ 
The  findings  reported  here  also  complicate  the  goa  o 
reaching  a  precise  atomic  description  of  rPA  protein  iso¬ 
form  content  and  its  relationship  to  rPA  product  stability 
by  dispelling  earlier  notions  of  a  proposed  direct  corre¬ 
spondence  between  site  specific  deamidation  and  isoform 
identity.  At  present,  we  showed  that  total  molecular  dea¬ 
midation  plainly  correlated  with  total  net  negative 
charge  and  isoform  complexity,  but  the  precise  combina¬ 
tions  of  site-specific  deamidation  among  the  different 
Asn  residues  and  visible  isoforms  remains  undefined.  In 
contrast,  the  loss  of  rPA  biological  activity  was  observed 
to  correlate  with  increased  deamidation  as  shown  next. 


Analysis  of  Biochemical  Activities:  Asn 
Deamidation  Associates  With  Loss  of 
Furin  Cleavage  and  Holotoxin  Assembly 

Three  interrelated  biochemical  activities  of  rPA  were 
measured  using  gel-shift  assays  to  investigate  a  possible 
effect  by  deamidation  on  protein  activity.  As  reported 
previously,21  these  assays  test  the  combined  functions  ot 
site-specific  proteolysis  by  furin  to  generate  PAes,  mter- 
molecular  protein  associations  that  create  the  heptamenc 
(PAe3)7  form,  and  then  binding  of  LF  to  create  the 
LFCPAeah  holoenzyme.  Modification  across  all  seven  Asn 
sites  at  a  mean  total  deamidation  ranging  from  10.5% 
(List  rPA)  to  18.2%  (rPA  lot  a,  1996  pool  A)  did  not  visi- 
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Fiq.  9.  Asparagine  deamidation  associates  with  reduced  biochemi¬ 
cal  activities  of  PA.  (A)  gel  shift  assay  for  furin-cleavage  activity.  (B)  ge 
shift  assay  for  holoenzyme-formation  activity.  Lot  identities  mean  dea¬ 
midation  of  rPA  (in  parentheses)  and  sample  content  per  lane  as  foi¬ 
ls  1  Ust  Biologicals  PA  (10.5%);  2.  GMP  rPA  <1Z2%fc 3  1 1996  pool 
A  rPA  (18  2%)’  4  1996  pool  B  rPA  (26.5%),  5,  1995  rPA  (27.2 /o),  , 

aged  GMP  rPA  (26.5%);  7° GMP  rPA;  8,  LF;  9,  LF  +  GMP  rP A;  10 ,  LF 
+91996  pool  A  rPA;  11,  LF  +1996  pool  B  rPA;  12,  LF  +  aged  GMP  rPA. 
Both  assays  performed  as  described  previously. 


formation  [Fig.  9(A  and  B),  lanes  2,  3,  9,  10,  and  legend]. 
On  the  other  hand,  full  inhibitions  of  these  activities 
occured  at  26.5%  mean  total  deamidation  [Fig.  9(A  and 
B),  lanes  4  and  11].  While  artificially  aged  GMP  protein 
containing  26%  mean  total  deamidation  retained  partial 
cleavage  activity,  it  did  not  retain  activity  for  holo¬ 
enzyme  formation  (Fig.  9,  lane  12).  Although  residues 
N713,  N537,  and  N719  had  the  greatest  incremental 
increase  in  deamidation  between  fresh  and  aged  protein 
[Fig.  10{A),  %  change],  the  apparent  correlation  cannot 
be  statistically  evaluated  and  therefore  we  cannot  yet 
propose  whether  deamidation  at  these  sites  exerts  a  pri¬ 
mary  causative  impact  on  loss  of  biochemical  function. 

Thus,  we  found  that  the  biochemical  and  biophysical 
activities  of  PA  were  sensitive  to  Asn  deamidation.  The 
ability  of  rPA  to  function  as  substrate  for  site-specific 
cleavage,  to  participate  in  heptamerization,  and  then 
bind  LF  to  form  holoenzyme  complex  were  inhibited 
above  a  threshold  level  of  26%  mean  total  deamidation 
across  the  seven  described  residues.  Although  deamida¬ 
tion  varied  among  individual  sites,  and  increased  differ¬ 
entially  by  artificial  aging,  their  pattern  of  deamidation 
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bly  affect  furin  cleavage,  heptamerization,  or  holoenzyme 
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. f'?.  I0'  c°roParis°n  of  computational  and  observed  deamidation  in  rPA.  (A)  Table  comparing  rank  order  of  observed  Asn  deamidation  to  that  pre- 
n  I  nlc  “ITT dea“on°r9  If  n9  database  identifier  for  PA  (1ACC).  Abbreviations  as  follows:  Asn  site,  amino  acid  coordinate 

°  t~?Uhre  D  1/2  h  f  J  ln  ,ys,°  deamida,ion  calculated  from  the  theoretical  coefficient  of  deamidation;  rank  order  for  propensity 

of  deamidation  as  listed  by.  f  predicted  rate  calculated  by  the  Robinson  algorithm;  F,  observed  percentage  for  fresh  PA  protein;  T,  observed  percent- 
nf9PA04trea,ed  Vh  n  L  i  n°  deamidation  observed;  (B)  Projection  of  seven  observed  sites  of  deamidation  onto  the  three-dimensional  structure 
of  PA.  Arrows  mark  deamidated  Asp  with  numeric  designation  by  coordinate  number,  as  defined  in  Figure  5;  pink  dots,  projection  of  original  aspar- 
gmyl  side  group.  Color  scheme  as  follows:  domain  1  (orange,  residues  1-258)  contains  two  calcium-binding  sites,  the  cleavage  site  foMurin  dass 
blndlf  for  EF  and  lethal  LF.  Domain  2  (red,  259-487)  forms  most  of  the  channel  and  protein-protein  interactions.  Domain  3 
p  ^  P  ay  V  °ll9°menzat'°n.  Domain  4  (blue,  596-735)  is  responsible  for  binding  the  cell-surface  receptor-  (C-  following 

afnfin  a  !  ,STh  observedA,sn  deamldation  sites  onto  a  space-filling  model  of  the  X-ray  structure  of  monomeric  PA,  and  rotated  180® 

nnm  .ho  J^9i^'S,  !  Vlewing.  Color  scheme  as  for  panel  B  except  domain  1  is  pink.;  (D)  Projection  of  seven  observed  Asn  deamidation  sites 

onto  the  structure  of  PA  heptameric  prepore  (PA^L,  central  pore  axis  angled  slightly  down  and  left,  and  one  LF-binding  domain  (pink) 

rerr’ov®d  (dotted  traPez0|d) ,  fo  expose  buried  sites  and  position  others  for  viewing.  Deamidation  at  residue  N537  is  highest  among  all  samples  tested 
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does  not  yet  allow  the  assignment  of  a  direct  causal  rela¬ 
tionship  between  the  loss  of  biochemical  activity  and 
modification  at  any  single  site.  Site-directed  mutagenesis 
to  change  individual  sites  to  Asp  may  help  to  elucidate  a 
specific  cause  for  loss  of  structural  activity. 

Analysis  of  Biological  Activity:  Cytotoxic  Activity 
Correlates  Inversely  With  Isoform  Complexity 
and  Asn  Deamidation 

We  also  employed  LFCrPA^L  holoenzyme,  formed  in 
the  same  manner,  to  test  cytotoxicity  as  an  indicator  of 
the  effect  of  deamidation  on  rPA  biological  activity.  This 
assay  is  believed  to  test  the  ability  of  PA  to  translocate 
LF  across  the  cell  membrane  which  then  kills  macro¬ 
phage  cells.  The  dose  of  protein  causing  50%  killing 


(ED50)  was  measured,  as  standardized  across  experi¬ 
ments  by  comparison  to  fully  active  PA.  For  GMP  grade 
rPA  and  R  &  D  lots  of  rPA,  both  ED50  and  site-specific 
deamidation  generally  increased  with  increasing  isoform 
complexity  (Figs.  5  and  6).  Comparing  mean  total  deami¬ 
dation  across  all  seven  measured  sites  to  ED50  values 
showed  a  statistical  correlation  between  deamidation 
and  loss  of  biological  activity  (r  =  0.96144,  P  =  0.0386) 
considering  the  four  lots  shown  in  Figure  6.  Statistically 
significant  correlation  between  the  loss  of  biological  ac¬ 
tivity  and  changes  in  individual  Asn  deamidation  were 
observed  as  follows:  N408,  r  =  0.99473,  P  =  0.0053; 
N466,  r  =  0.97894,  P  =  0.0211;  N601,  r  =  0.95546,  P  = 
0.0445.  Interestingly,  deamidation  in  the  most  reactive 
position,  N537,  did  not  vary  statistically  with  loss  of  bio¬ 
logical  activity  (r  =  0.88351,  P  =  0.1165).  Also,  although 
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Figure  10.  (Continued.) 


deamidation  at  the  second  most  deamidated  residue, 
N713,  appeared  to  vary  directly  with  loss  of  biological  ac¬ 
tivity,  statistical  parameters  (r  =  0.94365,  P  =  0.0563) 
fell  just  outside  levels  of  significance  and  therefore  a  de¬ 
finitive  association  was  not  supported.  Of  the  remaining 
sites  lacking  statistical  relationship  in  this  regard,  one 
showed  moderate  deamidation  (N719)  and  the  other 
showed  the  lowest  measured  level  of  deamidation  (N602) 
and  did  not  change.  Among  these  data,  only  deamidation 
levels  at  N408,  N466,  and  N601  varied  directly  and  sig¬ 
nificantly  with  the  loss  of  rPA  cytotoxic  activity. 

Thus,  we  found  no  supportive  evidence  of  deamidation 
in  rPA  protein  at  position  N162  which  juxtaposes  the 
furin  cleavage  site,63  but  these  data  indicated  that  inhi¬ 
bition  of  cleavage  correlated  with  deamidation  else¬ 
where.  In  spite  of  the  high  deamidation  at  N537,  and 
accelerated  deamidation  at  N713  and  N719  caused  by 
freeze/thaw/hold  treatment,  variation  in  none  of  these 
residues  correlated  statistically  with  loss  of  biological  ac¬ 
tivity.  Instead,  the  loss  of  activity  appeared  to  be  associ¬ 
ated  with  one  or  more  deamidations  at  positions  N408, 
N466,  and  N601  since  these  did  vary  with  statistical 
significance  upon  loss  of  biological  activity.  It  is  reasona¬ 
ble  that  the  sites  of  slowest  observed  degradation 


exerted  the  greatest  effect  on  biological  activity  and  vice 
versa,  otherwise  selective  pressure  would  have  gener¬ 
ated  different  structure/function  relationships  to  main¬ 
tain  anthrax  holotoxin  that  degrades  naturally  while 
retaining  biological  function  in  vivo.  This  hypothesis 
may  be  testable  using  site-directed  substitutions  in 
whole  PA  protein. 

With  regard  to  eliciting  protective  immunity  as  a  vac¬ 
cine  antigen,  we  note  that  although  both  of  the  major 
isoforms  provide  equivalent  protection  against  lethal  an¬ 
thrax  challenge  as  compared  to  the  whole  pharmaceuti¬ 
cal  preparation,23'70  we  now  know  that  neither  of  these 
primary  isoforms  represented  native  (unmodified)  pro¬ 
tein:  both  contained  various  amounts  of  deamidation  at 
the  seven  labile  positions  described  herein.  These  find¬ 
ings  demonstrate  that  certain  specific  deamidations 
present  in  the  test  article  did  not  negatively  impact  its 
measurable  efficacy  as  a  vaccine  substance  or  in  PA- 
mediated  biological  activity  in  vitro.  While  the  state  of 
protein  modification  in  vivo  during  disease  has  not  been 
reported,  PA  is  likely  to  degrade  at  these  residues  based 
on  our  findings  that  rPA  deamidates  rapidly  in  B. 
anthracis  during  expression.  Thus,  the  biological  activity 
of  rPA  is  apparently  tolerant  to  a  range  of  Asn  deamida- 
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final  site  (N537)  is  situated  in  domain  3  [Fig.  10(A)]. 
Residues  N601  and  N602  are  at  the  beginning  of  domain 
4,  buried  within  the  heptamer  molecule,  and  contact 
domains  2  and  3  of  the  same  PA  monomer  subunit.  Resi¬ 
dues  713  and  N719  are  exposed  on  outer  side  of  domain 
4  and  flank  the  ATR  binding  site  at  714-719.  A  struc¬ 
tural  rationale  for  why  deamidation  at  N601  correlates 
strongly  with  loss  of  biological  activity  was  not  immedi¬ 
ately  apparent.  Residue  N537,  the  most  highly  deami¬ 
dated  Asn  of  PA,  appears  also  to  be  the  most  exposed 
and  exterior  Asn  in  the  modeled  structures  [Fig.  10(D)], 
Its  positioning  at  the  apex  of  a  solvent-exposed  outer 
loop  [Fig.  10(A)]  may  also  support  its  high  tendency  to 
deamidate.  The  peripheral  location  of  this  residue  in  the 
heptamer,  away  from  ATR  or  effector  binding  regions, 
also  concurs  with  the  absence  of  a  prominent  effect  by 
N537  deamidation  on  the  biochemical,  biological,  and 
immunological  activity  of  rPA. 


CONCLUSIONS 


ture  or  stability.  Isoform  complexity  and  asparagine  dea¬ 
midation  associates  inversely  with  PA  biochemical  and 
biological  activity  at  four  residues:  N408,  N466,  N601, 
and  N713.  Deamidation  at  N408  and  N466  show  the 
strongest  statistical  correlation  with  the  loss  of  PA  bio¬ 
logical  activity  in  vitro.  Thus,  with  respect  to  measuring 
PA  degradation  by  deamidation,  N537  is  a  suitable  bio¬ 
marker  for  general  structural  stability,  while  N408  and 
N466  are  better  suited  for  assays  of  biological  stability. 
The  high  level  of  deamidation  at  N537  in  all  isoforms 
and  lack  of  apparent  effect  on  biochemical  or  biological 
activity  concurs  with  findings  from  our  prior  study 
which  demonstrated  that  the  two  primary  isoforms  of 
rPA  are  equally  protective  as  vaccine  immunogens 
against  lethal  anthrax  exposure.23,70  Finally,  the  method 
described  herein  for  collecting  and  analyzing  spectral 
data  is  generally  applicable  for  the  study  of  deamidation 
in  other  proteins  because  it  employs  an  LC-ESI-MS/MS 
instrument  of  nominal  mass  accuracy  and  because  calcu¬ 
lations  are  performed  strictly  from  direct  empirical  meas¬ 
urements  without  dependence  on  theoretical  models. 


We  draw  the  following  conclusions  from  this  study  of 
structural  heterogeneity  in  rPA.  Micro-charge  heteroge¬ 
neity  of  the  PA  protein  is  caused  by  partial  and  dissimilar 
deamidation  across  at  least  seven  asparagine  residues: 
N408,  N466,  N537,  N601,  N602,  N713,  and  N719.  For 
freshly  thawed  pharmaceutical  grade  rPA,  these  deami¬ 
dations  are  the  only  observed  post-translational  modifi¬ 
cations  capable  of  imparting  additional  charge  to  the 
protein,  and  no  other  amino  acid  modification  is  present 
except  for  a  minute  amount  of  methonine  oxidation 
(M521)  Position  N537  has  the  highest  cumulative  level 
of  deamidation  in  GMP  grade  rPA  protein,  commercial 
PA  protein,  research  and  development  grade  protein,  as 
well  as  in  purified  rPA  isoforms,  and  remains  the  high¬ 
est-ranking  modified  position  after  intentional  degrada¬ 
tion  by  freeze/thaw/hold  treatment.  Neither  of  the  two 
predominant  isoforms  present  in  GMP-grade  and  com¬ 
mercial  grade  PA  has  a  molecular  charge  corresponding 
to  native  protein,  and  both  of  these  isoforms  are  par¬ 
tially  deamidated  across  all  seven  of  the  labile  Asn.  No 
specific  site  of  deamidation  can  yet  be  attributed  to  the 
appearance  of  a  specific  isoform  band,  spot,  peak,  or  pat¬ 
tern.  By  inspection  of  isoform  patterning,  and  by  in¬ 
ference  from  measured  modifications  on  fresh  and  aged 
protein,  rPA  appears  to  be  rapidly  deamidated  m  vivo 
and  then  deamidates  to  a  smaller  incremental  degree  in 
vitro.  Mean  percent  deamidation  increases  in  direct  pro¬ 
portion  with  isoform  complexity.  After  intentional  degra¬ 
dation  by  freeze/thaw/hold  treatment,  the  rank  order  of 
percent  deamidation  among  these  seven  residues  does 
not  change,  but  the  apparent  rates  of  deamidation  for 
N713  and  N719  trends  toward  their  predicted  rates. 
Whether  measured  in  freshly  thawed  or  treated  protein, 
most  of  these  deamidations  do  not  concur  with  their 
theoretical  ordinal  rank  of  deamidation.  The  Robinson 
algorithm  for  predicting  deamidation  at  Asn  and  Gin 
residues  in  protein  is  presently  not  well  suited  for  prese¬ 
lecting  specific  residues  with  which  to  assay  rPA  struc- 
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